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Abstract 

This  report  presents  the  final  results  of  the  project  “An  improved  model  for  operational 
specification  of  the  electron  density  structure  up  to  geosynchronous  heights”  (FA5209-09-P- 
0253)  with  the  following  main  objectives: 

a)  To  upgrade  the  TaD  model 

b)  To  validate  the  results  using  ISR  and  topside  electron  density  profiles,  GNSS  TEC  and 
partial  TEC  parameters 

c)  To  specify  added-value  products  based  on  the  TaD  functional  capabilities 

d)  To  implement  the  TaD  model  on-line  ingesting  autoscaled  ionospheric  parameters  from 
Athens  Digisonde 

e)  To  develop  a  prototype  web  service  demonstrating  the  functionality  of  TaD  to  operate  on¬ 
line  using  as  input  autoscaled  ionospheric  parameters  and 

f)  To  demonstrate  some  of  the  indicative  value-added  products  derived  from  the  TaD 
outputs. 

Introduction 

Description  of  the  problem 

The  safety  and  security  of  space  operations  requires  that  spacecraft  operators  know  and 
understand  the  environment  around  their  spacecraft.  A  measure  of  the  ionospheric  element  of 
that  environment  is  given  by  the  precise  values  of  the  total  electron  content  (TEC)  between 
the  spacecraft  and  some  location  on  the  surface  of  the  Earth  (partial  TEC)  -  and  in  some  cases 
between  two  spacecrafts.  The  partial  TEC  value  differs  from  those  derived  from  ground-based 
GNSS  receivers.  The  latter  values  represent  the  TEC  of  the  whole  ionosphere-plasmasphere 
system,  i.e.  height-integrated  from  the  surface  to  the  GNSS  satellites  around  20000  km 
altitude.  For  many  space  activities  we  must  partition  GNSS  TEC  above  and  below  the 
spacecraft  in  order  to  derive  the  TEC  applicable  to  a  particular  activity. 

As  direct  observations  of  electron  density  are  sparse,  precise  estimation  of  the  partial  TEC 
requires  better  modelling  of  the  topside  ionosphere,  i.e.  the  region  above  peak  electron  density 
and  that  gradually  merges  into  the  plasmasphere  at  altitudes  above  the  O  -  H+  transition 
height.  The  topside  is  the  region  that  contains  most  of  the  plasma  that  contributes  to  the  TEC 
values  that  we  seek  to  estimate;  typically  only  25%  is  below  the  F2  peak  while  typically  < 
10%  comes  from  the  plasmasphere.  Furthermore,  the  topside  is  the  region  in  which  most  LEO 
spacecraft  operate,  so  accurate  estimates  of  the  TEC  relevant  to  a  particular  spacecraft  activity 
require  us  to  partition  the  TEC  within  the  topside.  Thus  modelling  the  topside  electron  density 
is  central  to  our  key  objective. 

Along  this  direction,  Kutiev  and  Marinov  (2007)  and  Kutiev  et  al.  (2006  and  2007)  developed 
an  empirical  model  of  the  O  -  H+  transition  height  (hi),  the  topside  electron  density  scale 
height  (fix)  and  their  ratio  Rt=HT/hT,  named  Topside  Sounder  Model  (TSM)  based  on  the 
Alouette/ISIS  database.  To  further  approach  the  problem  of  developing  a 
topside/plasmaspheric  electron  density  reconstruction  model  with  increased  accuracy,  Kutiev 
et  al.  (2009a)  offered  analytical  formulas  for  obtaining  the  shape  of  the  vertical  plasma 
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distribution  in  the  topside  ionosphere  and  plasmasphere  based  on  TSM  parameters.  This 
profiler  was  named  Topside  Sounder  Model  Profiler  (TSMP),  which  models  separately  the  0+ 
and  H+  density  profiles.  To  obtain  the  density  distribution,  it  needs  specification  of  the  F  layer 
maximum  density  (NmF2),  its  height  (hmF2)  and  its  scale  height  (FIm)  at  its  lower  boundary. 
Recently,  Kutiev  et  ah,  (2009b)  applied  TSMP  to  Digisonde  measurements  to  reconstruct  the 
electron  density  profile  from  the  F  layer  peak  to  GNSS  orbits  (TaD  model).  A  new  expression 
of  FI+ scale  height  in  the  plasmasphere  was  extracted  from  ISIS-1  topside  sounder  data,  as  a 
function  of  geomagnetic  latitude.  The  verification  of  the  TaD  model  has  been  attempted  by 
Belehaki  et  al.  (2009a),  comparing  the  integral  of  TaD  profiles  with  total  electron  content 
(TEC)  measured  by  the  CHAMP  satellite  (from  400  km  up  to  GPS  orbits)  and  by  ground- 
based  GPS  receivers.  The  fitting  between  the  measured  GPS-TEC  and  the  TaD-TEC  depends 
a  lot  from  the  geographic  location  of  the  test  site  while  the  comparison  with  CHAMP  TEC 
data  showed  some  deviations  that  could  be  also  attributed  to  the  background  models  applied 
to  reconstruct  the  CHAMP  electron  density  profiles.  Overall,  the  results  from  the  first 
verification  analysis  revealed  the  need  for  further  improvements,  which  is  the  main  objective 
of  this  project. 


Summary  of  the  main  results 

In  the  frames  of  this  project  the  Topside  Sounders  Model  assisted  by  Digisonde  (TaD)  is 
further  improved  (Section  1)  according  to  the  following  methodology:  a)  Improvement  of  the 
TaD  scaling  technique  based  on  the  calculation  of  0+,  H+,  and  He+  density  distributions  in 
transition  region  between  topside  F  region  and  plasmasphere,  extracted  from  the  analysis  of 
the  electron  density  profiles  from  ISIS-1.  This  yield  a  more  reliable  determination  of 
important  scale  height  parameters  (0+,  H+,  and  He+)  and  transition  height;  b)  Improvements 
of  the  TSMP  formulation  to  achieve  optimization  of  the  TaD  algorithm.  The  main  model 
expressions  were  revised  to  include  He+  distribution  as  a  function  of  geomagnetic  latitude  and 
local  time. 

A  systematic  validation  of  the  new  profiler  is  presented  in  Section  2  using  Topside  and  ISR 
Electron  Density  Profiles  (EDP),  TEC  and  partial  TEC  parameters.  The  results  show 
reduction  of  the  model  error  2.25  times  in  comparison  to  the  previous  version  of  the  TaD 
model.  Validation  of  the  model  results  based  on  comparison  with  ISIS  1  EDP  shows  clearly 
the  model’s  ability  to  reproduce  with  impressive  accuracy  the  ISIS  1  EDP  (98.8%). 
Systematical  comparison  between  the  0+  distribution  of  TaD  and  measured  EDP  from  the 
Malvern  ISR  gave  a  mean  RMSE  of  12%  while  this  error  increased  slightly  when  compared 
between  the  total  electron  density  and  the  observed  values.  A  discussion  related  to  these 
results  is  given  in  the  last  part  of  this  Section. 

Based  on  TaD  capabilities  and  on  users’  requests  we  have  been  able  to  specify  and  develop  a 
set  of  value  added  products,  available  through  this  prototype  web  service,  while  additional 
products  have  been  specified  (Section  3)  for  future  implementation  in  the  DIAS  system. 

The  results  from  the  optimized  TaD  algorithm  are  given  in  Section  4,  and  are  available  on¬ 
line  through  the  web  address  http  ://www.iono  .noa.  gr/ElectronDensity/EDProfile  .php .  For  this 
on-line  implementation,  autoscaled  data  from  Athens  Digisonde  are  ingested  in  the  algorithm, 
to  simulate  real-time  operation  conditions.  A  technical  failure  in  Athens  Digisonde  prevented 
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us  to  run  in  real-time  this  prototype  web  service.  The  Athens  Digisonde  is  currently  out  of 
operation  due  to  technical  failure  caused  by  a  severe  thunderstorm  with  lightning  activity  that 
hit  Penteli  area  on  Saturday  3  July  2010.  The  functionality  of  the  prototype  in  real-time 
conditions  will  be  resumed  as  soon  as  the  Athens  Digisonde  becomes  fully  operational  again. 
Currently  this  prototype  web  service  provides  data  and  value  added  products  for  four  years 
period:  22  July  2006  to  2  July  2010. 

In  our  attempt  to  validate  the  TaD  results  and  demonstrate  the  functionality  of  the  algorithms 
in  real-time,  we  have  specified  a  number  of  issues  that  should  be  considered  in  the  future, 
especially  when  this  model  is  going  to  be  integrated  to  the  DIAS  system,  and  develop 
products  based  on  real-time  data  from  a  number  of  ionospheric  stations,  such  as  maps  of  TEC. 
These  issues  among  others  are  discussed  in  Section  5:  Future  steps. 
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1.  Development  of  the  advanced  model  for  the  electron  density 
reconstruction  up  to  geosynchronous  height 

1.1  Background  models 

The  development  of  the  advanced  model  for  the  electron  density  reconstruction  up  to 
geosynchronous  height  accomplished  within  this  project  continued  a  series  of  successive 
approximations  to  the  open  issue  of  the  reconstruction  of  the  electron  density  profile  up  to 
plasmaspheric  heights.  The  improved  technique  is  based  on:  i)  the  Topside  Sounders  Model 
(TSM),  ii)  the  Topside  Sounder  Model  Profiler  (TSMP)  and  iii)  the  TSMP-assisted  Digisonde 
(TaD)  profile. 

The  Topside  Sounders  Model  (TSM)  reproduces  the  0+-  H+  (upper)  transition  height  (hT) 
and  the  topside  electron  density  scale  height  (HT),  based  on  172,622  measured  Ne  profiles  by 
topside  sounders  on  Alouette-la,  -lb,  -lc  and  -2  and  ISIS-1  and  -2  satellites  (Bilitza,  2001). 
The  Ht  and  hT  models  are  represented  by  5 -dimensional  polynomials  expressing  both 
quantities  as  functions  of  month  of  the  year,  geomagnetic  latitude,  local  time,  solar  flux  F107, 
and  the  geomagnetic  index  Kp.  Later,  Kutiev  and  Marinov  (2007)  developed  a  model  of  the 
ratio  Rt=Hx/hT,  based  on  ratios  of  individual  profiles  and  combined  the  three  models  into  a 
single  model  named  Topside  Sounder  Model  (TSM).  The  key  element  in  the  modeling  is  the 
assumption  that  the  lowest  gradient  on  the  topside  Ne  profile  represents  the  altitude  gradient 
of  the  O  density  profile.  Assuming  exponential  distribution  of  0+  density,  the  transition 
height  hT  is  found  at  the  height,  where  the  upward  extrapolated  0+  density  becomes  one  half 
of  the  measured  Ne  density.  The  0+  gradient,  inferred  from  the  Ne  profile  is  converted  into 
scale  height,  e.g.  the  distance  in  km  at  which  the  0+  density  changes  e  (Euler’s  number) 
times.  Therefore,  both  Hj  and  hi  have  the  same  dimension.  TSM  provides  the  topside  F 
region  scale  height  H  i  and  transition  height  hi  between  0+  and  H+  dominant  regions,  based  on 
measured  data  only.  As  these  parameters  are  intrinsic  to  plasma  distribution  in  ionosphere  and 
plasmasphere,  they  can  be  used  in  any  model  dealing  with  the  plasma  density  structure.  It  has 
to  be  noted,  that  Ht  is  a  vertical  scale  height  and  differ  from  the  theoretical  plasma  scale 
height  defined  along  magnetic  field  lines  (see  Belehaki  et  ah,  2006). 

Validation  tests  of  TSM  performance  include  the  comparison  with  two  of  the  most  popular 
ionospheric  profile  reconstruction  techniques,  the  Parametirized  Ionospheric  Model  -  PIM 
(Daniell  et  ah,  1995)  and  the  NeQuick  model  (Hochegger  et  ah,  2000;  Radicella  and 
Leitinger,  2001).  The  results  provide  usefull  input  of  the  model’s  performance  regarding  the 
season  and  the  solar  activity  level  (Stankov  et  al.,  2007). 

Topside  Sounder  Model  Profiler  (TSMP):  Most  recently,  Kutiev  et  al.  (2006)  offered  a 
method  for  obtaining  the  shape  of  the  vertical  plasma  distribution  in  topside  ionosphere  and 
plasmasphere  by  introducing  Ht  and  hx  parameters  into  three  well  known  formulas  describing 
the  vertical  plasma  distribution:  a-Chapman,  Sech-squared,  and  Exponential  (see  Stankov  et 
al.,  2003).  The  analytical  formulas  containing  HT  and  hT  were  named  as  Topside  Sounder 
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Model  Profiler  (TSMP).  TSMP  provides  the  shape  of  the  vertical  plasma  distribution;  the 
density  distribution  is  obtained  once  the  F  layer  density  (NmF2)  and  height  hmF2  at  its  lower 
boundary  is  specified. 

TSMP  provides  the  shape  of  the  topside  N(h)  profile  as  the  sum  of  the  0+  and  H+  profiles 
from  hmF2  up  to  GPS  heights.  The  0+  profile  takes  HT  as  its  own  scale  height  (Ho+=  HT),  at 
transition  height  hT  both  densities  equalize,  and  above  H+  scale  height  (Hh+)  is  taken  as  16 
times  0+  scale  height.  Later,  Kutiev  et  al.  (2009b;  2009c)  found  that  the  Hh+/Ho+  (or  Hp/HT) 
ratio  is  not  constant  with  a  value  of  16,  but  depends  on  latitude  and  on  average  ranges 
between  12  and  4.  The  present  formulation  of  TSMP  has  the  form: 


f 

Ne  =  N(),  (h)  +  N(y  (hT  )  exp 

v 


h  -  hT  | 
Hp  , 


(1) 


where: 


N0+{h)  =  Nm  exp< 


h~hm 


f 

exp 

v 


h-hP 

H  + 

o*  ) 

Here  Ne  is  the  electron  density  assumed  to  be  equal  to  the  sum  of  0+  and  H+  densities  (Ne  = 
N0++  Nh+)  in  all  altitudes.  N0+  is  described  by  a-Chapman  function,  with  a  scale  height  H0+ 
being  one  half  of  Ht,  while  Nh+  is  pure  exponential  function  with  a  scale  height  Hh  h  denoted 
as  Hp.  Below  hT,  NH+  is  forced  to  decrease  by  taking  absolute  difference  |h-  hT  |  in  order  to 
comply  with  the  real  distribution.  Therefore,  the  maximum  of  Nh+  distribution  is  fixed  at 
transition  height,  although  in  reality  it  lies  slightly  above.  At  transition  height,  Nh+  =  No+(  h  r) 
and  scale  the  hydrogen  ion  distribution  in  plasmasphere.  To  calculate  the  topside  density 
profile,  TSMP  needs  values  of  NmF2  and  hmF2  at  the  lower  boundary  (denoted  as  Nm  and 
hm). 

TSMP-assisted  Digisonde  (TaD)  profile:  As  mentioned  above,  TSM  provides  also  the  ratio 
Rt=HT/hT,  which  is  not  directly  used  in  TSMP.  The  model  ratio  Rt,  which  actually  over¬ 
determine  the  Ne  profile,  possess  an  important  property:  it  can  be  used  as  a  link  between  TSM 
parameters  and  another  measured  topside  quantities.  This  property  is  used  by  Kutiev  et  al. 
(2009a)  to  link  TSMP  to  Digisonde  topside  profiling  technique.  The  new  hybrid  profiler, 
named  TSMP-assisted  Digisonde  (TaD)  profiler,  uses  ionogram-derived  scale  height  at  hmF2 
(denoted  as  Hm),  corrected  in  a  way  to  comply  with  TSM-produced  HT.  This  correction  factor 
depends  on  the  location,  but  in  conducted  analysis  an  average  value  of  1.25  was  taken. 
Therefore  the  corrected  scale  height  is  calculated  as  Hd  =  1.25Hm.  Further,  TaD  uses  the 
model  ratio  Rt  to  obtain  its  actual  transition  height  hTD  multiplying  Hd  by  Rt.  It  was  found 
that  the  plasmasphere  scale  height  Hh+  (now  denoted  as  Hp)  is  not  equal  to  16  times  Ht,  as 
assumed  from  theoretical  considerations,  but  it  is  latitude  dependent,  ranging  in  average  from 
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about  12  at  equator  to  5  at  80°  geomagnetic  latitude.  An  empirical  relation  was  introduced, 
defining  HP  =  (9.cos2(glat)+4)Hd,  where  glut  was  geomagnetic  latitude. 

In  TaD  approach,  the  lower  boundary  values  and  the  scale  height  are  provided  by  Digisonde 
measurements.  In  other  words,  TaD  feeds  TSMP  with  Digisonde  data.  This  combination 
makes  TaD  profiling  connected  with  measurements  and  offers  variety  of  applications, 
especially  the  reconstruction  of  actual  Ne  profiles  over  each  Digisonde  location.  This  property 
of  TaD  can  be  directly  implemented  in  any  regional  digisonde  network  such  as  the  European 
DIAS  system  to  provide  3-D  distribution  (latitude,  longitude,  altitude)  of  electron  density  in 
real  time 

Belehaki  et  al.  (2009a)  conducted  an  extensive  work  on  verification  of  TaD  profiling 
technique  by  comparing  the  integral  of  TaD  profiles  with  total  electron  content  (TEC) 
measured  by  CHAMP  satellite  (from  400  km  up  to  GPS  orbits)  and  ground-based  GPS 
receivers  at  locations  of  Athens  and  Juliusruh.  Comparison  of  TaD  with  CHAMP-derived 
TEC  showed  standard  deviation  of  4.6  and  3.4  TECU  and  correlation  coefficient  of  0.86  and 
0.81  for  Athens  and  Juliusruh  respectively.  Comparison  with  ground-based  TEC  yielded  a 
standard  deviation  of  5.3  and  2.9  TECU  and  correlation  coefficient  of  0.92  and  0.95  for 
corresponding  locations. 


1.2.  Development  of  the  improved  TaD  model 
1.2.1.  Improvement  of  the  scaling  technique 

In  the  Interim  report  (Belehaki  et  al.,  2009b),  we  stressed  attention  on  the  problem  of  how 
well  TaD  profiler  describes  the  transition  region  between  the  topside  F  region  and 
plasmasphere.  We  found  that  a  certain  number  of  Ne  profiles  measured  by  ISIS-1  topside 
sounder  require  additional  amount  of  light  ions  to  that  of  H+  in  order  to  equalize  0+  density  at 
transition  height  h  r.  We  assumed  that  these  Ne  profiles  indicated  the  presence  of  He+  ions 
(Heelis  et  al.,  1990;  Su  et  al.,  2005)  and  showed  the  latitude  and  local  time  distribution  of  the 
He+  density,  shown  in  Figure  1.  We  also  included  the  helium  ion  density  in  the  TSMP 
formula  describing  the  electron  density  as  a  function  of  altitude  by  introducing  a  new 
parameter  g  being  the  ratio  of  H+  and  0+  densities,  g  =  n(H+)/n(0+).  TSMP  formula  is  now 
given  in  the  form: 


f 

Ne  =  N()l  ( h )  +  gN()  (h,  ) exp 

v 


h  -  hT  ' 
Hp  , 


f 

+  (1  ~g)N0+  {hT )  exp 

v 


\h-hT[ 
4  Ht  , 


(2) 


In  the  absence  of  He+,  g  =  1  and  the  third  term  is  equal  to  zero.  The  lower  is  g,  the  more  He+ 
ions  are  present.  Obviously,  g  should  be  a  function  of  geomagnetic  latitude,  local  time  and 
probably  to  the  season.  It  is  reasonable  to  assume  that  He+  has  its  maximum  again  at  the 
transition  height,  hj,  and  exhibits  a  pure  exponential  distribution  above,  with  a  scale  height 
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equal  to  4  times  that  of  Hr.  However  at  higher  latitudes,  for  the  plasmaspheric  scale  height  Hp, 
which  theoretically  should  be  16  times  that  of  Hx,  this  factor  drops  significantly  due  to  the 
shape  of  plasmasphere.  Having  in  mind  this  fact,  here  we  try  to  infer  the  factor  s  =  Hhc+/Ht 
from  the  data.  We  further  denote  He+  scale  height  (HHe+)  as  He.  So,  formula  (2)  is  now 
modified  as  in  last  exponent  instead  of  4HT  we  place  sHT  and  the  new  TSMP  formula  is  given 
by  the  following  equation: 


f 

Ne  =  N()l  ( h )  +  gN0+  (hT ) exp 

V 


h  —  hT  | 
Hp  , 


f 

+  (l-#ot(^r)ex  P 

V 


\h-hT[ 

sHt 


(3) 
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Figure  1:  The  geomagnetic  latitude/local  time  contour  map  of  the  difference  of  0+  and  H+ 
densities  at  transition  height. 


By  definition  He+  density  strongly  depends  on  O  and  H+  densities  at  transition  height,  which 
in  their  turn  strongly  depend  on  scale  heights  Ht  and  Hp,  detennined  from  lowest  and  highest 
altitude  gradients  of  the  measured  Ne  profiles.  Increasing  the  accuracy  of  fitting  the  measured 
Ne  profiles  in  transition  region  comes  up  through  the  increase  of  the  accuracy  of  He+ 
determination  at  hT.  To  increase  the  accuracy  of  Hp  and  HT,  we  developed  a  sophisticated 
algorithm,  based  on  optimization  of  HT,  Hp,  and  s  towards  best  fitting  with  measured  Ne 
profile.  Formula  (3)  is  run  multiple  times  by  using  a  set  of  HT,  Hp,  and  5  and  the  model 
profile  at  each  run  is  compared  with  measured  Ne  profile  by  calculating  the  root  mean  square 
(RMS)  error.  A  series  of  values  is  assigned  for  each  of  these  variables  and  Ne  is  calculated 
from  formula  (3),  executing  three  nested  cycles:  the  outer  for  Hx,  middle  for  Hp  and  internal 
for  s.  Each  run  takes  one  of  the  H  r  values,  calculates  the  corresponding  hx,  takes  one  of  the 
Hp  values,  calculates  the  ratio  g  at  hT,  takes  one  of  the  5  values,  and  then  calculates  based  on 
formula  (3)  the  whole  profile  in  the  altitude  range  corresponding  to  the  measured  profile.  The 
RMS  error  is  stores  and  the  next  cycle  is  performed.  We  have  set  40  values  around  the 
minimum  Hx,  40  values  around  maximum  Hp,  and  120  values  of  s  around  4.  The  whole 
optimization  of  a  given  profile  contains  40  x  40  x  120  =  192,000  runs.  So,  we  select  the  set  of 


Page  11  of  44 


AOARD  Contract  No  FA5209-09-P-0253,  Final  Report 


the  three  variables  that  has  the  smallest  RMS  error.  In  this  way  we  perform  optimization  to 
14,775  measured  ISIS-1  profiles.  Figure  2  shows,  as  a  sample,  the  profile  measured  on  day 
97  of  year  1969,  17:03  UT  around  equator  (lat  =  12°S,  long=12.5°E).  Logarithm  of  ion 
density  profiles  obtained  without  optimization  (denoted  as  “old”)  and  those  with  optimization 
(denoted  by  “new”),  are  color  coded,  along  with  the  measured  Ne  profile  presented  with  the 
black  line.  The  “old”  profiles  are  presented  by  dashed  lines,  while  the  “new”  profiles  are 
given  by  solid  lines.  The  total  ion  density  presented  the  sum  of  the  “new”  profiles  is  given  by 
the  orange  line.  In  this  sample,  optimization  increases  slightly  the  values  of  Hj  and  hj  and 
reduces  Hp.  Optimization  invokes  its  largest  effect  on  He+  density,  which  is  substantially 
reduced.  However,  He+  density  profile  below  transition  height  hT  depends  strongly  on  H+ 
density  shape,  which  as  we  have  demonstrated,  could  not  represent  the  real  situation.  As  we 
are  interested  in  reconstructing  the  electron  density  profile,  which  below  hj  is  mainly  defined 
by  0+,  the  shape  of  H+  is  not  important.  The  same  is  true  for  the  He+  profile  below  hT.  In  this 
study  we  are  mainly  interested  in  He+  density  around  its  maximum  fixed  at  hT.  The  standard 
deviation  of  the  optimized  model  from  measured  profiles,  denoted  as  RMS,  is  0.051,  while 
that  without  optimization  is  0.112  units  of  logarithm  density.  Here  optimization  has  reduced 
RMS  error  more  than  twice. 

We  found  that  the  He+  scale  height  cannot  be  reliably  obtained.  Indeed,  the  RMS  error  is 
formed  by  comparing  formula  (3)  with  the  whole  measured  profile.  He+  density  in  some  cases 
is  comparable  with  electron  density  within  the  transition  region,  but  it  is  always  much  smaller 
outside  that  region  and  its  net  contribution  to  the  RMS  error  is  negligible.  When  calculating 
procedure  executes  formula  (3)  for  variables  in  the  nested  cycles,  RMS  errors  are  practically 
determined  by  the  pairs  of  HT  and  Hp  values,  with  no  significance  of  values  of  5.  Thus  5  value 
could  not  be  optimized.  Statistically  for  the  whole  database,  s  values  are  spread  almost 
uniformly  in  the  defined  range.  Taking  into  account  this  fact,  we  further  take  s  =  4. 
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Figure  2:  A  sample  profile  measured  on  day  97  of  year  1969,  17:03  UT  around  equator  (lat  = 
12°S,  long=12.5°E).  Logarithm  of  ion  density  profiles  obtained  without  optimization  (denoted 
as  “old”)  and  those  with  optimization  (denoted  by  “new”),  are  color  coded,  along  with  the 
measured  Ne  profile  presented  with  the  black  line.  The  “old”  profiles  are  presented  by  dashed 
lines,  while  the  “new”  profiles  are  given  by  solid  lines.  The  total  ion  density  presented  the 
sum  of  the  “new”  profiles  is  given  by  the  orange  line. 

The  histogram  of  distribution  of  RMS  errors  collected  from  the  whole  selected  ISIS-1 
database  containing  14,775  profiles  is  given  in  Figure  3.  The  pink  filled  bars  represent 
histogram  of  RMS  without  optimization,  while  transparent  blue  bars  show  distribution  of 
RMS  error  after  optimization.  The  average  RMS  error  without  optimization  is  0.189,  while 
the  average  RMS  error  with  optimization  is  0.085.  Optimization  reduced  in  average  the 
model  error  2.25  times. 
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Figure  3:  A  histogram  of  distribution  of  RMS  errors  collected  from  the  whole  selected  ISIS-1 
database  containing  14,775  profiles.  The  pink  filled  bars  represent  histogram  of  RMS  without 
optimization,  while  transparent  blue  bars  show  distribution  of  RMS  error  after  optimization. 


1.2.2  Improvements  of  the  TSMP  formulation 

The  g  function  determines  the  presence  of  He+  at  transition  height  hi  and  it  is  a  key  parameter 
for  the  reconstruction  of  the  electron  density.  TaD  cannot  obtain  any  information  about  the 
He+  from  the  bottomside  sounding  and  therefore  the  presence  of  helium  ions  is  a  priori 
determined  on  the  base  of  topside  sounder  profiles.  For  obtaining  analytical  expression  for  g, 
we  reduced  the  ISIS-1  database,  excluding  profiles  with  g  =  1.  The  latter  restriction  is 
necessary  in  order  to  exclude  contribution  of  the  profiles  with  no  presence  of  He+.  To 
approximate  g  obtained  from  individual  measured  profiles  (g  =  n(H+)/  n(0+)),  we  chose 
analytical  function  of  three  variables:  local  time,  geomagnetic  latitude  and  O  density.  The 
latter  variable  was  added  to  the  first  two  after  specifying  the  strong  correlation  between  g  and 
n(0+).  The  inclusion  of  n(0+)  in  analytical  expression  significantly  reduced  the  model  error, 
e.g.  standard  deviation  between  model  and  measured  values. 

Mathematical  formulation 

Coefficient  matrix  CFN  is  a  solution  of  the  problem  for  LSQ-approximation  of  the  given  data 
points  (. x(k);f(k )),  x{k)={xi{k),  xi{k),  xfk),),  i.e.  CFN  minimizes  the  functional: 
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fj(/(k)-F(CFN,x(k)))2 

k= 1 


(4) 


where  F(CFN,x )  =  F(x)  is  an  element  of  the  LSQ-approximation. 

We  define  the  function  F(xijC2-F3)  within  the  intervals  [a„h(],  \c„di\,  i=  1, _ ,3;  base-functions 

and  coefficients  as  follows: 

F(x i ,  x2 ,  x3 )  =  J  R, .  (.s’, )  £  fi2/.  (.s2  (.s3  )CFN(j\ ,  j 2 ,  y3 ) 

ti=i  A=1  h=  1  PI 


where: 

(xpXj^JeR3  ,  (s[,s,,s3)eR3  fori  =  \,2,2> 

xt  G  [a,  A]  >  ^  G  s,  =C  +T^(<  -c/)» 

b~a, 

CFN{ix ,  i2 , /, )  are  elements  of  coefficient  matrix  with  a  size  nl  x  n2  x  «3 ; 

Base-functions  are  Bj  . (5,. ),  j  =  1, 2,...,ni . 

For  the  base-function  we  use: 

Algebraic  functions:  1,5,52,..., s*,...; 

Tchebishev’s  functions:  To(j)=1,  Ti(j)=j,...,  Tk(.s)=2.sTk-i(s)  -  Tk-2(s),  for  k=2,3, 
i.e.  Tk(^)=cos(k  arccos(.s)); 

Trigonometric  functions:  1,  sin(.s),  cos(.s),  sin(Cs),  cos (ks), 

In  our  case  the  base-function  and  coefficient,  defined  above,  are: 

"l  =7;  [«i,&iJ=[o,24];  [c1,d1J=[o,2^];  (s)=1 ,  >2(s')=sin(s),51 3(.?)=cos(s)  (.s)=cos(3.s) 

«2=7;  [«2,&2]=  [-90,90];  [c2,d2]=[-y,y];  B2  l(s)  =  l ,  B2 j2(s)  =  sin(s),  ...,B2 ?  (.v)=sin(6.s). 

»3  =5;  [a3,Z>3]=[4,12];  [c3,c?3]=[- l,l] ;  B^ ^(s)=l  ,B3 j2(s)=j,  B^  3(s)=2.?2 -1  ;^3  ^(s)  =  4.v3  -  3s... 
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Figure  4:  The  scatter  plot  of  model  g  ( gm )  versus  measured  g  values.  The  thick  blue  line  is 
the  linear  regression  through  the  origin  over  gm  values. 

The  scatter  plot  of  model  g  (gm)  versus  measured  g  values  is  given  in  Figure  4.  The  thick 
blue  line  is  the  linear  regression  through  the  origin  over  gm  values.  Correlation  coefficient  is 
0.98,  which  shows  that  the  analytical  function  of  g  represents  the  data  exceptionally  well. 
Function  g  is  a  3-D  surface  that  could  not  be  easily  visualized.  To  get  an  idea  how  g  function 
depends  on  each  of  the  three  variables,  in  Figure  5  we  have  approximated  the  projection  of 
gm  values  on  each  of  axes  plane  with  low  order  polynomials.  In  the  top  panel,  the  linear  fit 
shows  that  gm  have  a  tendency  to  decrease  with  increasing  0+  density.  This  means  that  He+ 
predominantly  appears  (reversely  proportional  to  g)  when  O H  density  increase,  which  happens 
at  the  F  region  anomaly  crests.  The  middle  panel  shows  that  gm  has  a  minimum  of  around  0.6 
at  low  latitudes  and  increases  toward  higher  latitudes,  which  complies  with  the  known  fact 
that  He+  is  predominantly  present  around  equator.  Local  time  dependence  shows  a  minimum 
in  the  early  afternoon  hours,  when  the  equatorial  anomaly  is  well  developed.  As  a  conclusion, 
we  may  say  that  gm  behavior  is  in  compliance  with  expected  variation  of  He+  density  and  the 
function  can  reliably  estimate  contribution  of  helium  ions  in  transition  region. 


Page  16  of  44 


AOARD  Contract  No  FA5209-09-P-0253,  Final  Report 


0  3  6  9  12  15  18  21  24 

local  time 

Figure  5:  The  projection  of  the  function  gm  on  each  of  axes  plane  with  low  order 
polynomials;  Top  panel:  the  linear  fit  shows  that  gm  have  a  tendency  to  decrease  with 
increasing  0+  density;  Middle  panel:  gm  has  a  minimum  of  around  0.6  at  low  latitudes  and 
increases  toward  higher  latitudes;  Bottom  panel:  local  time  dependence  shows  a  minimum  in 
the  early  afternoon  hours. 
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2.  Validation  of  the  new  profiler 


Validation  is  a  strong  requirement  for  implementation  of  a  model  for  operational  use.  Model 
validation  is  the  comprehensive,  systematic  quantitative  comparison  of  model  output  with 
observations  and  is  required  for  identifying  and  documenting  model  strengths  and  weaknesses. 
For  the  validation  of  the  TaD  profiler  we  have  performed  a  systematic  comparison  of  the 
model  results  with  all  available  ground  and  space  based  observations.  These  datasets  were 
identified  through  the  state-of-art  survey  carried  out  within  the  framework  of  this  project 
(Belehaki  et  ah,  2009b)  and  include  TEC  data  from  ground-based  GNSS  receivers,  partial 
TEC  data  from  CHAMP  mission  and  ISR  EDP. 

2.1  Comparison  with  measured  profiles  from  ISIS-1 


As  a  first  test  of  the  performance  of  the  new  profiler,  we  compared  TaD  electron  density 
profiles  with  actual  measurements  obtained  from  ISIS-1  satellite.  This  is  not  an  independent 
evaluation  test  as  ISIS-1  profiles  were  used  for  the  development  of  the  new  method,  but 
provides  with  evidence  of  how  well  the  new  profiler  reproduces  the  input  measurements.  The 
comparisons  are  made  through  the  estimation  of  the  RMSE  over  whole  profile,  according  to 
the  following  formula: 


where  up  indicates  the  number  of  points  in  each  ISIS  1  EDP,  hmin  and  hmax  are  the  lowest  and 
highest  altitudes  for  which  ISIS  1  electron  density  measurements  are  provided  for  a  given 
profile,  nr  ad  and  n/sisi  are  the  electron  density  values  at  a  given  height  h,  provided  by  the 
TAD  model  and  ISIS  1  satellite  respectively. 
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Figure  6:  The  histogram  of  distribution  of  RMSE  between  ISIS  1  measured  electron  density 
and  TaD  output  values. 
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The  histogram  of  distribution  of  RMSE  is  presented  in  Figure  6,  indicating  a  very  successful 
fit  between  measurements  and  model  predictions.  This  very  good  agreement  can  also  be  seen 
between  the  TEC  parameters  estimated  by  the  ISIS  1  profiles  and  those  provided  by  the  TaD 
model.  In  Figure  7a  we  present  the  correlation  between  the  TEC  parameter  estimated  for  the 
plasmasphere  region,  integrating  the  electron  density  from  the  transition  height  up  to  the 
maximum  height  of  the  ISIS  1  satellite.  A  very  high  correlation  coefficient  of  98.2%  is 
estimated  between  the  measured  and  the  modeled  parameters.  Similarly,  in  Figure  7b  we 
present  the  TEC  parameter  for  the  topside  ionosphere,  integrating  from  the  lower  height  of  the 
ISIS  1  satellite  (slightly  above  hmF2)  up  to  the  transition  height.  Here  the  correlation 
coefficient,  between  observed  and  modeled  quantities,  is  also  impressively  high  reaching  the 
99%. 


Figure  7:  Scatter  plots  between  the  ISIS  1  TEC  parameter  and  the  TaD  TEC  for  (a)  the 
plasmasphere  region,  integrating  the  electron  density  from  the  transition  height  up  to  the 
maximum  height  of  the  ISIS  1  satellite  and  (b)  the  topside  ionosphere,  integrating  from  the 
lower  height  of  the  ISIS  satellite  (slightly  above  hmF2)  up  to  the  transition  height. 
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2.2  Comparison  with  TEC  derived  from  ground-based  GNSS  receivers 

In  this  section,  the  comparison  is  based  on  TEC  parameters  estimated  from  data  broadcasted 
by  GNSS  receiving  stations  and  correspond  to  Athens  (38.000N,  23.500E)  and  Juliusruh 
(54.600N,  13.40E)  coordinates.  The  method  of  comparison  is  the  following:  the  bottomside 
integral  of  Ne  profile  obtained  from  Digisonde  ionograms  is  added  to  the  integral  of  topside 
Ne  profile  obtained  by  TaD  method  and  the  sum  is  compared  with  TEC  values  obtained  by 
GPS  signals  at  the  same  locations.  The  database  is  compiled  from  ionosonde  and 
corresponding  GPS-TEC  data  at  30  min  sampling  compassing  5  years  (2001-2005)  period.  To 
give  the  degree  of  agreement  between  the  two  products,  the  GPS-derived  TEC  and  the  TaD- 
derived  TEC,  we  present  their  scatter  plot  for  Athens  and  Juliusruh  sites,  in  Figures  8  and  9 
respectively.  With  green  we  mark  the  regression  line. 


Athens  2001-2005 


120  - 


0  40  80  120 

TEC  (GPS)  in  TECU 


Figure  8:  Comparison  between  GPS-derived  TEC  and  TaD-derived  TEC  over  Athens. 
Regression  line  is  also  plotted  in  green. 


Page  20  of  44 


AOARD  Contract  No  FA5209-09-P-0253,  Final  Report 


0  40  80  120 
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Figure  9:  Comparison  between  GPS-derived  TEC  and  TaD-derived  TEC  over  Juliusruh. 
Regression  line  is  also  plotted  in  green. 


It  is  important  here  to  give  a  measure  of  the  improvement  of  the  TaD  after  its  optimization. 
Table  1  gives  the  average  value  of  the  difference  dTEC: 

dTEC  =  TEC  (TaD)  -  TEC(GPS) 

on  the  whole  sample  of  TEC  parameters  for  each  location,  with  the  standard  deviation  of  the 
population. 


Tablel:  Improvement  of  the  optimized  TaD  performance  in  respect  to  the  original  version, 
estimated  for  Athens  and  Juliusruh  sites  for  the  period  2001-2005 


Athens 

Juliusruh 

Original  TaD 

Optimized 

TaD 

Original  TaD 

Optimized 

TaD 

No  of  points 

63,317 

71,088 

dTEC  mean  value 
(TECU) 

-1.9 

-1.5 

-1.0 

-0.8 

dTEC  Standard 
Deviation  (TECU) 

5.3 

3.4 

2.9 

1.9 
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2.3  Comparison  with  data  from  Malvern  ISR 

An  important  part  of  the  improved  TaD  validation  plan  is  the  comparison  of  the  model 
profiles  with  measured  profiles  obtained  from  Malvern  ISR.  This  test  provides  infonnation  for 
the  model’s  performance  in  the  altitude  zone  200-600  km  (approximately).  The  importance  of 
this  test  is  threefold:  i)  ISR  profiles  are  the  only  source  of  full  electron  density  observation 
from  the  ground,  ii)  they  consist  a  totally  independent  source  of  infonnation  and  iii)  this 
comparison  is  strongly  suggested  by  US  National  Space  Weather  Program  and  recent  studies 
for  the  establishment  of  the  top  high  priority  metric  for  the  ionosphere-thermosphere  domain. 
For  this  purpose,  a  data  set  of  some  of  4000  electron  density  profiles  obtained  at  Malvern  ISR 
site  (52.1  °N,  2.3  °E)  during  the  time  interval  1968-1971  were  analyzed.  After  quality  check 
almost  the  30%  of  profiles  have  been  eliminated  from  our  sample.  The  F  layer  peak  density 
NmF2  and  height  hmF2  were  extracted  from  each  individual  profile.  However,  for  the  scaling 
of  the  ISR  EDPs  we  have  applied  the  method  used  in  the  topside  sounder  profiles  scaling,  i.e. 
the  NmF2  and  hmF2  parameters  have  been  located  in  each  measured  profile;  then  its  topside 
scale  height  HTm  has  been  obtained  as  a  regression  line  from  a  number  of  data  points  above 
the  peak,  with  the  transition  height  hTM  being  approximated  as  the  height  where  0+  is  one  half 
of  the  measured  density.  The  four  quantities  taken  from  the  measured  profile  have  been 
ingested  to  the  formula  (3)  and  the  model  profile  was  calculated  for  altitude  range  of  the 
topside  part  of  measured  profile.  Keeping  this  in  mind,  the  analysis  attempted  below  may  be 
considered  as  a  rough  evaluation  of  the  method’s  perfonnance  but  in  our  opinion  it  provides 
suggestive  conclusions. 

Figure  10  gives  some  examples  of  the  derived  profiles. 

Malvem-ISR  profiles  cover  the  altitude  range  up  to  approx.  700  km  which  is  about  the 
altitude  of  the  actual  transition  height,  so  it  is  expected  that  0+  will  dominates  and  the 
observed  Ne  profile  will  follow  closely  the  TaD  profile  of  the  0+  density.  Based  on  this,  we 
first  compared  the  model  predictions  for  0+  distribution  in  respect  to  the  measured  profiles. 
The  results  are  expressed  in  terms  of  the  nonnalized  RMSE  estimates  that  where  obtained 
over  the  topside  profiles  (from  hmF2  up  to  the  maximum  height  of  the  observed  profiles). 
Their  distribution  is  presented  in  Figure  11. 

We  also  compared  the  observed  Ne  profiles  with  models  predictions  for  the  total  density.  The 
corresponding  results  for  the  normalized  RMSE  are  presented  in  Figure  12. 
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In(density) 


In(density) 


Figure  10:  Examples  for  TaD  derived  profiles  based  at  Malvern  site.  The  ISR  EDP  are 
denoted  with  the  cross  symbol.  Red  line  indicates  the  modeled  0+  profile  and  blue  line 
indicates  the  modeled  EDP. 
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Normalised  RMSE  (%):  Topside  profile  (0+) 
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Figure  11:  The  normalized  RMSE  distribution  of  model  predictions  for  0+  in  respect  to  the 

measured  profiles. 
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Figure  12:  The  normalized  RMSE  distribution  of  modeled  total  electron  density  in  respect  to 

the  measured  profiles. 


The  comparison  between  Figure  11  and  12  confirms  our  assumption  for  better  fitting  of  the 
TaD  0+  profiles  with  the  measured  EDP. 

Due  to  the  fact  that  the  topside  ISR  profiles  are  limited  to  700  km,  transition  height  cannot  be 
always  reliably  obtained,  e.g.  we  cannot  reach  the  height  where  0+  is  one  half  of  the  measured 
density.  Moreover,  the  measured  profiles  around  their  upper  values  exhibit  an  artificial 
increase  which  TaD  wrongly  is  interpreted  as  appearance  of  H+. 

This  might  be  the  reason  of  the  tendency  of  the  model  to  overestimate  the  electron  density,  as 
seen  in  Figure  13.  Here  we  present  the  distribution  of  the  simple  difference  between  the 
observed  and  the  modeled  electron  densities  at  three  selected  altitudes:  400  km  (top),  500  km 
(middle)  and  600  km  (bottom).  Another  interesting  result  concerning  the  dependence  of  the 
model’s  performance  with  the  height  is  that  the  fit  of  the  model  to  the  observed  values 
increases  with  increased  altitude. 
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Observed  -  Modelled  (400  km) 


Observed  -  Modelled  (500  km) 


Observed  -  Modelled  (600  km) 


Figure  13:  The  distribution  of  the  simple  difference  between  the  observed  ISR  and  the 
modeled  electron  densities  at  400  km  (top),  500  km  (middle)  and  600  km  (bottom). 
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As  TEC  is  the  most  critical  parameter  for  the  reliable  performance  of  several  applications,  we 
also  performed  comparison  tests  between  TEC  estimates  obtained  from  modeled  and  observed 
EDP  from  the  ISR  at  Malvern  site  concentrating  in  the  topside  part  of  the  profile. 

First,  we  present  in  Figure  14  the  scatter  plot  of  TEC  modeled  versus  TEC  observed 
estimates.  The  results  indicate  that  the  two  parameters  correlate  reasonably  well.  A  larger 
scatter  is  seen  in  the  daytime  hours  (the  area  with  larger  TEC  values),  which  is  probably  due 
to  the  high  altitude  of  the  transition  height,  something  that  makes  its  accurate  determination 
problematic  due  to  the  lack  of  observations  above  700  km.  The  O  /H+  transition  height  varies 
but  seldom  drops  below  500  km  at  night  or  800  km  in  the  daytime,  although  it  may  lie  as  high 
as  1100  km,  depending  on  the  geophysical  conditions,  and  particularly  on  solar  activity 
(Denton  etal.,  1999). 


Malvern,  1968-1971 


Figure  14:  The  scatter  plot  of  TEC  modeled  versus  TEC  observed  estimates 

In  Figure  15  we  present  the  TEC  error  (abs(TECobs-TECmod))  distribution  and  the  relative 
TEC  error  distribution.  From  the  top  to  the  bottom  we  present  the  results  based  on  all 
available  measurement,  in  the  middle  we  present  the  corresponding  distribution  diagrams  for 
a  subset  of  measurements  that  correspond  to  daytime  EDP  and  in  the  bottom  we  give  the 
distribution  diagrams  for  nighttime  cases.  The  analysis  on  the  whole  sample  of  measurements 
gives  s  mean  error  of  3TECU  and  mean  relative  error  of  28.5%.  However  this  error  reduces 
considerably  when  we  include  only  the  nighttime  cases.  This  is  an  additional  evidence 
showing  that  during  daytime,  when  the  transition  height  is  in  its  upper  limits,  or  even  higher 
of  700  km,  there  are  not  enough  data  for  the  TaD  to  estimate  correctly  the  transition  height 
and  this  yields  to  artificial  results.  The  same  problem  was  reported  by  Belehaki  and  Kersley 
(2006),  where  they  compared  the  ISR  EDP  from  Malvern  with  the  results  from  the  topisde 
extrapolation  based  on  the  Reinisch  and  Huang  model  (2001). 
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Distribution  of  TEC  error  (in  TECU)  Distribution  of  TEC  relative  error  (%) 
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Figure  15:  The  TEC  error  (abs(TECobs-TECmod))  distribution  and  the  relative  TEC  error 
distribution.  From  the  top  to  the  bottom:  a)  results  based  on  all  available  measurements;  b)  the 
corresponding  distribution  diagrams  for  a  subset  of  measurements  that  correspond  to  daytime 
EDP;  c)  the  distribution  diagrams  for  nighttime  cases. 
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2.4  Summary  on  validation  results 


•  The  TaD  model  reproduces  with  very  high  accuracy  the  ISIS  1  profiles  demonstrating  that 
the  model  does  what  it  is  designed  to  do  and  that  it  is  free  of  intrinsic  uncertainties. 

•  The  model  error  of  3TECU  is  close  to  the  measurement  (GNSS)  error 

•  The  mean  RMSE  obtained  from  the  comparison  of  ISR  EDPs  is  in  general  larger  than  the 
RMSE  obtained  from  the  comparison  with  ISIS-1  profiles  but  this  was  expected  as  ISIS-1 
profiles  were  also  used  for  the  model’s  development.  The  comparison  with  ISR  EDPs 
gives  a  model’s  error  mainly  attributed  to  limitations  imposed  by  the  measured  profiles 
themselves.  During  nighttime  hours,  when  the  transition  height  is  lower  than  700  km  and 
can  be  detennined  by  the  TaD  algorithm,  the  model  values  are  in  very  high  correlation 
with  the  ISR  observations,  with  a  mean  error  of  only  1.14  TECU. 
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3.  Definition  of  value-added  products  based  on  TaD  model 


The  specification  of  value-added  products  is  based  on  TaD  functional  capabilities,  on  the 
survey  of  related  products  available  by  various  organizations  worldwide  and  on  latest  needs 
specified  by  various  users  groups,  identified  in  the  frames  of  the  European  COST  Action 
ES0803  “Developing  Space  Weather  Products  and  Services  in  Europe”. 

The  functional  capabilities  of  the  TaD  include:  a)  real-time  calculation  of  the  electron  density 
up  to  GEO  heights,  at  locations  where  a  digisonde  operate  b)  real-time  calculation  of  TEC, 
and  of  its  ionospheric  and  plasmaspheric  portions  c)  real-time  calculation  of  slab  thickness 
and  of  its  ionospheric  part  d)  real-time  availability  of  the  partial  TEC,  whose  altitude  ranges 
may  be  defined  by  the  user. 

In  addition  the  survey  of  related  products  indicated  that  worldwide  there  are  several 
organizations  providing  TEC  related  products,  however,  we  could  not  locate  any  web  site  of 
open  access  providing  data  of  partial  TEC.  More  precisely,  the  Space  Weather  Application 
Center  in  DLR,  Germany,  receives  GPS  data  and  calculates  the  following  European  products: 
TEC  Map,  Error  Map,  Rate  of  Change,  Longitudinal  Gradient,  Latitudinal  Gradient.  The 
SPECTRE  service  in  France,  based  on  the  same  GPS  data  provides  TEC  maps,  relative  error 
maps  and  electronic  biases.  They  also  provide  related  STEC  and  VTEC  products.  In  Japan, 
three  research  centers  (Kyoto  University,  Nagoya  University  and  the  National  Institute  of 
Information  and  Communications  Technology)  provide  TEC  maps  over  Japan,  over  North 
America  and  over  Europe  based  on  GPS  data.  Global  maps  of  TEC  are  provided  by  the 
MIT/Haystack  Observatory  in  USA  and  by  NASA  in  the  standard  IONEX  format. 

Latest  needs  specified  by  various  users  groups  indicated  that  for  a  number  of  space  activities 
we  must  partition  GNSS  TEC  above  and  below  the  spacecraft  in  order  to  derive  the  TEC 
applicable  to  a  particular  activity.  Such  activities  are  the  radar  tracking  of  spacecraft  as  part  of 
space  surveillance  programmes;  the  use  of  space-based  synthetic  aperture  radars  (SAR)  to 
carry  out  remote  sensing  of  the  Earth’s  surface;  the  use  of  on-board  GNSS  for  small  satellite 
missions. 

Based  on  the  above  analysis  it  is  evident  that  a  model  such  as  TaD  able  to  provide  the 
analytical  function  of  the  electron  density  up  to  GEO  height  would  support  a  number  of 
operational  applications.  As  the  operation  of  this  model  relies  on  the  availability  of 
ionospheric  characteristics  of  the  F-region,  in  Europe  the  TaD  model  could  be  integrated  to 
the  DIAS  system  (the  European  Digital  Upper  Atmosphere  Server,  http://dias.space.noa.gr)  in 
order  to  make  available  a  comprehensive  set  of  services  characterizing  conditions  in  the 
topside  ionosphere  and  the  plasmasphere.  However  that  set  up  of  such  a  comprehensive 
service  requires  several  steps  including  a)  application  of  the  TaD  to  as  many  as  possible 
Digisonde  locations  b)  development  of  mapping  techniques  over  an  area  (such  as  Europe)  c) 
computation  of  the  3D  grid  of  ED  values  over  the  specified  area. 

Within  the  limited  time  of  this  project,  we  have  been  able  to  demonstrate  the  operation  of  TaD 
using  as  test  site  Athens,  where  a  Digisonde  operate  and  provides  in  real-time  autoscaled 
ionospheric  characteristics. 
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The  flow  chart  of  the  set  of  web  services  providing  the  TaD  output  and  the  related  value 
added  products  is  shown  in  Figure  16. 
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Figure  16:  Flow  chart  of  TaD  web  services 
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All  TaD  data  and  products  are  calculated  with  15  min  resolution.  The  list  of  all  available 

TaD  related  web  services  is  given  below: 

•  The  analytical  function  of  the  electron  density  profile  up  to  20,000  km  in  ASCII  values 

•  Visualization  of  simultaneous  ionograms  and  their  corresponding  TaD  electron  density 
extrapolated  profile. 

•  The  TEC  up  to  20,000  km,  in  daily  plots 

•  The  partial  TEC  corresponding  to  the  plasmaspheric  heights  (from  the  transition  height  up 
to  20,000  km),  in  daily  plots 

•  The  partial  TEC  corresponding  to  the  topside  ionosphere  (from  the  hmF2  up  to  transition 
height),  in  daily  plots 

•  The  slab  thickness,  in  daily  plots 

•  The  ionospheric  part  of  the  slab  thickness,  in  daily  plots 
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4.  On-line  demonstration  of  the  TaD  results  and  value  added 
products 

A  prototype  service  has  been  set  up  to  demonstrate  the  functionality  of  the  TaD  algorithm  and 
to  distribute  some  indicative  value-added  products. 

The  service  is  available  through  the  address: 


http://www.iono.noa.gr/ElectronDensity/EDProlile.php 


The  web  service  has  been  developed  with  the  following  open  source  tools: 

a)  PHP,  for  the  web  services  set  up, 

b)  Perl,  for  the  on-line  implementation  of  TaD  the  algorithm  that  has  been  converted 
from  its  original  FORT  AN  version 

c)  MySQL  for  the  database,  which  was  based  on  the  enhancement  of  the  DIAS  database. 
Below  we  present  some  screen  shots  to  demonstrate  the  functionality  of  the  web  service. 
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Figure  17:  The  Home  Page  of  the  TaD  web  services.  Area  (1)  includes  the  main  menu  for  quick  search  and  for  the  download  of  the  TaD  EDP  in  ASCII.  Area 
(2)  includes  the  visualization  of  the  TaD  EDP  and  its  corresponding  ionograms.  Area  (3)  includes  the  navigation  buttons 
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Figure  21:  The  analytical  function  of  the  electron  density  profile  up  to  20,000  km  in  ASCII 
values  downloadable  from  the  main  menu  of  the  Home  Page 
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5.  Further  steps 


As  a  next  step  we  propose  the  implementation  of  TaD  in  DIAS  system 
(http://dias.space.noa.gr)  to  make  available  the  TaD  EDP  at  all  DIAS  locations  (Athens, 
Rome,  Ebro,  Arenosillo,  Pruhonice,  Chilton,  Juliusruh,  Moscow)  and  consequently  to  be  able 
to  provide  important  services,  such  as,  maps  of  TEC  and  of  its  time  derivative,  maps  of 
plasmaspheric  TEC,  maps  of  slab  thickness,  3D  maps  of  electron  density  over  Europe.  This 
plan  requires 

a)  application  of  the  TaD  to  as  many  as  possible  DIAS  locations  and  if  possible  to  additional 
sites 

b)  development  of  mapping  techniques  over  an  area  (such  as  Europe) 

c)  upgrade  of  the  DIAS  user  layer  profile 

For  the  operational  release  of  the  system  the  following  important  issues  must  be  considered 
related  to  the  quality  of  real-time  autoscaled  ionospheric  data: 

-  Automatic  scaling  errors  affecting  the  TaD  perfonnance 

-  Data  gaps 

-  Definition  of  the  reference  level  in  order  to  determine  a  disturbance  index 

The  implementation  of  this  plan  is  subject  to  the  availability  of  the  necessary  funding 
resources. 
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